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Vortex Imaging in the 7r-Band of Magnesium Diboride 
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We report scanning tunneling spectroscopy imaging of the vortex lattice in single crystalline MgB2 . 
By tunneling parallel to the c-axis, a single superconducting gap (A = 2.2 meV) associated with the 
7r-band is observed. The vortices in the 7r-band have a large core size compared to estimates based 
on H C 2, and show an absence of localized states in the core. Furthermore, superconductivity between 
the vortices is rapidly suppressed by an applied field. These results suggest that superconductivity 
in the 7r-band is, at least partially, induced by the intrinsically superconducting cr-band. 

PACS numbers: PACS numbers: 74.50.+r, 74.70.Ad,74.60.Ec 
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Superconductivity in magnesium diboride (MgB2) 
with a remarkably high T c = 39 K was recently reported 
by Nagamatsu et al. [0. Since then, great attention has 
been directed towards understanding the detailed nature 
of superconductivity in this material, an in particular 
whether this is a one- or two-gap superconductor. Two- 
gap superconductivity was predicted theoretically j|, [|, 
and is now supported by an increasing number of exper- 
imental reports §, || | || gj], (yj || . Two-gap or 
two-band superconductivity was hrst studied in the fifties 
p3| , and has now found renewed relevance in MgB2. In 
addition, and contrary to many materials or alloys stud- 
ied earlier, the two bands in MgB2 have roughly equal 
filling factor, opening the possibility for interesting new 
phenomena. However, the exact microscopic details are 
still largely unexplored. An ideal way to address this is- 
sue is by local spectroscopic investigations of the mixed 
state, which has become possible with the recent avail- 
ability of high quality MgB 2 single crystals. 

In this Letter we report on scanning tunneling spec- 
troscopy (STS) measurements on single crystal MgB2, 
including the first vortex imaging in this material. Tun- 
neling parallel to the c-axis, we are able to selectively 
measure only the 7r-band |jl4| , in which the vortices are 
found to have a number of remarkable properties: An 
absence of localized states, a very large vortex core size 
compared to the estimate based on -ff c2 , and a strong 
core overlap. 

The STS experiments were performed using a home 
built scanning tunneling microscope (STM) installed in 
a 3 He, ultra high vacuum cryostat holding a 14 T mag- 
net |l5| ]. The measurements were done on the surface of 
an as grown single crystal, using electrochemically etched 
iridium tips. Single crystals of MgB 2 were grown using a 
high pressure method described elsewhere [ fl6| , yielding 
platelike samples with the surface normal parallel to the 
crystalline c-axis. The surface of the crystals are roughly 
0.25 x 0.25 mm 2 , and the thickness of the order of mi- 
crons. The critical temperature is typically T c w 38 — 39 
K, with a sharp transition, AT C = 0.5 K, measured by 



SQUID magnetometry |17j|. The STS experiments were 
done with both the tunneling direction and the applied 
magnetic field parallel to the c-axis, and the differen- 
tial conductivity measured using a standard AC lock-in 
technique. In this configuration the upper critical field 
extrapolates to H c2 (T = OK) = 3.1 T Using the 

Ginzburg-Landau (GL) expression for H c2 = (j) /(2Tr£ i 2 ), 
where <pQ = h/2e is the flux quantum, yields a coher- 
ence length, £gl = 10 nm. An estimate of the mean free 
path, based on the measured residual resistivity jlq] and 
specific heat [11 , and the calculated Fermi velocity p4| , 
gives Z = 50 — 100 nm, indicating that the samples are in 
the clean limit. 

We will first focus on the zero-field electronic spec- 
trum of MgB 2 . The observation of a single or double 
gaps depends on the orientation of the sample, as shown 
by Iavarone et al., who investigated a number of single 
grains with different, but unknown absolute orientations 
|12f . Here we report the first STS measurements on a 
MgB 2 single crystal, which allow a correlation between 
the tunneling direction and the observed gap(s). In Fig. 
1 we show a superconducting spectrum obtained at a 
temperature of 320 mK. This is an average of 40 spec- 
tra obtained along a 100 nm path, which shows perfect 
homogeneity. One observes a single gap with coherence 
peaks at ±2.5 meV, and additional weak shoulders at 
±6 meV, as indicated by the arrows. In addition, the 
flat region around the Fermi energy proves the absence 
of nodes in the gap, and hence that MgB 2 is a s-wave 
superconductor. The very low zero bias conductance in- 
dicates a high quality tunnel junction and a low noise 
level. True vacuum tunneling conditions were assured by 
varying the tunnel resistance, i?t, and verifying that the 
spectra normalized to the conductance outside the super- 
conducting gap collapse on a single curve. The spectrum 
can be fitted by the BCS expression for the density of 
states (DOS), including a finite quasiparticle lifetime, T 
]l9[ , and an experimental broadening. The result of the 
fit is shown in Fig. 1 and yields a superconducting gap, 
A = 2.2 meV. We have studied the temperature depen- 
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FIG. 1: Zero-field superconducting spectrum of MgB2 at 320 
mK for tunneling parallel to the c-axis and a tunnel resistance, 
Rt = 0.2 GO. (U = 0.1 V; I = 0.5 nA). The bias voltage is 
applied to the sample. Clear coherence peaks are seen at ±2.5 
meV, and additional weak shoulders at ±6 meV as indicated 
by the arrows. The line is a fit to the Dynes DOS Q ( A = 2.2 
meV, r = 0.1 meV) convoluted by a gaussian of width 0.5 
meV RMS to account for experimental smearing, including 
the use of a finite AC excitation (0.4 mV RMS). 



dence of the superconducting gap and found excellent 
agreement with the BCS A(T). 

The fact that only one superconducting gap is observed 
for tunneling parallel to the c-axis, is explained theoreti- 
cally by calculations of the Fermi surface and an analysis 
of how tunneling along different directions is coupled to 
the different bands. The Fermi surface of MgB2 falls into 
two distinct sheets: One is derived from cr-antibonding 
states of the boron p xy orbitals and is a two-dimensional 
cylindrical sheet parallel to c* , while the other consists 
of 7r-bonding and antibonding states of the boron p z or- 
bitals and is three-dimensional EjL |2(|. The tunnel- 
ing matrix element is different for the two bands, and 
depends on the tunneling direction. The c-axis tunnel- 
ing probability into the cr-band is ten times smaller than 
into the 7r-band Q. Furthermore, the calculated super- 
conducting gap sizes for the two different Fermi surfaces 
are different, with A a 7 meV, and 2 meV |$). 

This is in agreement with our results, where one gap with 
A = Ajr = 2.2 meV is observed, with the shoulders at 6 
meV being a remnant of A CT . The selective sensitivity to 
A^ turns out to be particularly useful, as we will show 
in the following. 

We now turn to measurements in an applied magnetic 
field. In a type-II superconductor such as MgB 2 , a mag- 
netic field penetrates into the sample in the form of vor- 
tices each carrying one flux quantum, which are gener- 
ally arranged in a periodic array: the vortex lattice. In 
the core of each vortex, superconductivity is suppressed 
within a radius roughly given by the coherence length, £. 



The vortex spacing is determined by the applied field and 
the flux quantization, and in the case of a hexagonal vor- 
tex lattice it is d = (2/^3 cfto/H) 1 / 2 . The magnetic fields 
were applied at 2 K, and the system allowed to stabilize 
for at least a few hours. After this time no vortex motion 
was observed, indicating a fast relaxation and hence low 
vortex pinning in the crystal. 

In Fig. 2a we show a STS image of a single vortex 
induced by a field of 0.05 T. The image was obtained 
by measuring the differential conductance at zero bias 
and normalizing this to the conductance at the coher- 
ence peak. Low values of the normalized conductance 
correpond to superconducting areas, and high values to 
the vortex cores. The low field is equivalent to a sep- 
aration, d = 220 nm 3> £. The vortices can therefore 
be considered as isolated from each other. Such isolated 
vortices are expected to contain localized quasiparticle 
states, which should show up as a zero bias conductance 
(ZBC) peak at the vortex centre provided that the 
sample is sufficiently clean to prevent these to be smeared 
out by scattering. We have measured the evolution of 
the spectra at a large number of positions along a trace 
across the vortex core as shown in Fig. 2c. Contrary to 
expectations, we find that the normalized ZBC increases 
to one with no indication of any localized states. Instead, 
the spectra in the centre of the vortex are absolutely flat, 
with no excess spectral weight at or close to zero bias. 
This absence of localized states is striking, considering 
that I = 5 — 10 x £gl ■ However, as we will show be- 
low, the coherence length in the 7r-band is approximately 
50 nm. This is much larger that the estimate based on 
H C 2, and equal to only one to two times the mean free 
path. Nonetheless, systematic studies of Nbi_ 2; Ta a: Se2 
with x = — 0.2, showed that even for (/i « 1 some 
excess weight close to zero bias was observed p2]. In 
parallel to STS imaging, STM topographic images were 
recorded (not shown), which revealed a flat surface with 
a RMS roughness of 6 A over the whole image area. 

Before analyzing the vortex profile in detail, we will 
consider the situation at higher fields. In Fig. 2b we show 
the STS image of the hexagonal vortex lattice observed 
at 0.2 T. We notice that the normalized ZBC between 
the vortices, is now enhanced with respect to the value 
far from the single vortex at 0.05 T. This increase of the 
"bulk" ZBC is unusual at such a modest field, only about 
7% of H c2 . To elucidate this behaviour, "bulk" spectra 
for fields between zero and 0.5 T are shown in Fig. 2d. 
It is clear that even modest fields rapidly suppress su- 
perconductivity in the region between the vortices. This 
is seen both by an increase of the ZBC and by a sup- 
pression of the coherence peaks outside the vortex cores, 
which one would only expect for fields close to H C 2, corre- 
sponding to a significant core overlap. This is consistent 
with earlier point contact spectroscopy measurements [|J , 
with the addition that we resolve the local behaviour on 
a microscopic scale. 
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FIG. 2: (color) Vortices in MgB 2 . Top: 250 x 250 nm 2 false 
color spectroscopic images of a single vortex induced by an 
applied field of 0.05 T (a), and the vortex lattice at 0.2 T 
(b). In both cases the tunnel resistance was, Rt = 0.4 Gfl 
(U — 0.2 V; J = 0.5 nA). The conductance is normalized to 
respectively 2.9 meV (0.05 T) and 3.9 meV (0.2 T). (c): 250 
nm trace across the single vortex indicated by the white line 
in (a), with spectra recorded each 2 nm. Each spectrum is 
normalized to the conductivity at 12 meV, and Rt — 0.4 GQ. 
A spectrum at the vortex centre together with one far from 
the vortex core have been highlighted in red for clarity, (d): 
Normalized spectra measured in zero field, and between the 
vortices for fields between 0.05 T and 0.5 T (R t = 0.2 GQ). 
Each subsequent spectrum is offset by 0.5 with respect to the 
previous one. The bars at zero bias indicate the respective 
zero conductivity for the offset spectra. All measurements in 
this figure were performed at 2 K. 
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FIG. 3: Vortex size and core overlap, (a): Normalized zero 
bias conductance versus distance from the centre, for the iso- 
lated vortex shown in Fig. 2a. The line is a fit to eq. (|l|) 
in the text, (b): Calculated "bulk" ZBC (left axis) and elec- 
tronic specific heat, 7, for 'yZ/l'n = 0.55/0.45 (right axis). 
The calculated values are compared to respectively the mea- 
sured bulk ZBC (circles), and specific heat measurements on 
polycrystalline samples (squares) 



We now return to the single vortex measurement. In 
Fig. 3a, we have plotted the normalized ZBC, a'(x,0), 
for the vortex trace measured at 0.05 T. It is immediately 
clear that the spatial extension of the vortex core is much 
larger than the 10 nm estimated from H C 2- Somewhat 
surprisingly, we find that the ZBC profile can be fitted 
by one minus the GL expression for the superconducting 
order parameter: 



a'(x,0) 



+ (1 - a' Q ) x (1 - tanhx/O, 



(1) 



where <j' = 0.068 is the normalized ZBC measured in 
zero field. The fit, shown in Fig. 3a, yields a coherence 
length, £ = £n = 49.6 ± 0.9 nm. Using the GL expres- 
sion to calculate the upper critical field with this value 
of the coherence length, yields H' c2 = 0.13 T. At 0.2 T 
we hence find ourselves in the bizarre situation of imag- 
ing the vortex lattice above the nominal H' c2 ■ Additional 
vortex lattice imaging has been performed as high as 0.5 
T. 
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This apparent paradox can be reconciled if one assumes 
that superconductivity in the 7r-band is induced by the 
CT-band, either by interband scattering, or Cooper pair 
tunneling |l^, ^3). This means that isolated the 7r-band 
would either be non-superconducting or have a very low 
upper critical field. Consequently the observed behaviour 
reflects the state in the a-band by an interband proxim- 
ity effect, along the lines of recent theoretical work p3| . 
The vanishing of A n (T) at the bulk T c further supports 
this conclusion. Finally, it is also consistent with esti- 
mates of the coherence lengths, using the BCS expression 
£o = ^Up/(7tA(0)) and considering each band separately. 
Taking the calculated average Fermi velocity in the ab- 
plane for the 7r-band, Vp — 5.35 x 10 5 m/s fl4|] , and 
the measured gap value = 2.2 meV we get £J = 51 
nm, in excellent agreement with £ w obtained from the 
vortex profile. A similar analysis for the er-band, using 
v% = 4.4 x 10 5 m/s || and A^ = 7.1 meV |J) yields 
£o = 13 nm. This agrees with the coherence length ob- 
tained from H C 2, and reinforces the conclusion that it 
is mainly the er-band which is responsible for supercon- 
ductivity in MgB2, and thus determins the macroscopic 
parameters T c and -ff C 2- 

As described above, it is the transfer from the cr-band 
that makes superconductivity in the 7r-band possible, de- 
spite a strong vortex core overlap already at very low 
magnetic fields. Constructing a simple model for the core 
overlap in the 7r-band by 



a'(r, 0) = a' + (1 - a' Q ) X 1 - IT tanh 



, (2) 



where are the vortex positions for a hexagonal lattice 
with a density corresponding to the applied field, we can 
calculate the ZBC at any position in the vortex lattice 
unit cell. In Fig. 3b we compare the calculated conduc- 
tivity at the midpoint between three vortices with the 
measured bulk ZBC. This shows a very good agreement, 
especially taking into account that there are no free pa- 
rameters in the calculation: £„. is determined from the 
vortex profile, and a' from the zero field measurement. 
The vortex core overlap also explains the, strongly non- 
linear field dependence of the electronic specific heat, 7 
j|, [ll| . In strongly type-II superconductors core overlap 
is usually negligible, with each vortex creating the same 
number of quasiparticles at the Fermi surface, and hence 
contributing by the same amount to the specific heat. In 
that case 7 = 7 n H/H C 2, where 7 n is the electronic spe- 
cific heat in the normal state. However, in the case of 
MgB2, with strong core overlap in the 7r-band, the iso- 
lated vortex assumption is violated. Instead, one can cal- 
culate the contribution from the 7r-band, simply by aver- 
aging the normalized ZBC in one vortex lattice unit cell, 
7w = 7n ( cr '( r ;0)) p3[ . On the other hand, the a-band 
can be described by the usual linear field dependence 



7 ff = 7^ H/H c2 . Adding the terms gives 7 = 7^ + j a , 
where 7„/7„ * s the relative weigth of the two bands. The 
calculated field dependence of 7 is shown in Fig. 3b, in 
perfect agreement with the measured specific heat for 
polycrystalline MgB 2 0, using 7^/7^ = 0.55/0.45. 

In summary, we have presented STS data on the tt- 
band in MgB 2 , including the first vortex imaging in this 
material. We have demonstrated the absence of localized 
states in the vortex core, a very large vortex core size and 
a strong core overlap. The data presents a striking ex- 
perimental demonstration of the fundamentally different 
microscopic properties of the two bands in MgB2 . 
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